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Abstract 



Given a polynomial x G M" i-^ p{x) in n = 2 variables, a symbolic-numerical 
algorithm is first described for detecting whether the connected component of the 
plane sublevel set V = {x : p{x) > 0} containing the origin is rigidly convex, 
^ . or equivalently, whether it has a linear matrix inequality (LMI) representation, 

CN \ or equivalently, if polynomial p{x) is hyperbolic with respect to the origin. The 

' problem boils down to checking whether a univariate polynomial matrix is positive 

■ semidefinite, an optimization problem that can be solved with eigenvalue decompo- 

sition. When the variety C = {x : p{x) = 0} is an algebraic curve of genus zero, a 
CD I second algorithm based on Bezoutians is proposed to detect whether V has an LMI 

representation and to build such a representation from a rational parametrization 
of C. Finally, some extensions to positive genus curves and to the case n > 2 are 
mentioned. 
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1 Introduction 



Linear matrix inequalities (LMIs) are versatile modeling objects in the context of con- 
vex programming, with many engineering applications [5j. An n-dimensional LMI set is 
defined as 

n 

^ = {x G : F{x) = Fq + ^x^F^ h 0} (1) 
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where the Fi G R"^^™ are given symmetric matrices of size m and >z means positive 
semidefinite. From the characteristic polynomial 

t ^ det(t/„ + F{x)) = poix)+pi{x)t + ■ ■ ■ +p„_i(a;)t'^ + t"" 

it follows from e.g. [361 Theorem 20] that 

^ = {x G : po{x) > 0, . . .,Pm-i{x) > 0}. (2) 

Hence the LMI set is basic semialgebraic: it is the intersection of polynomial sublevel 
sets. From linearity of F{x) and convexity of the cone of positive semidefinite matrices, 
it also follows that JF is convex. Hence LMI sets are convex basic semialgebraic. 
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Figure 1: The TV screen level set is not LMI. 

One may then wonder whether all convex basic semialgebraic sets are LMI. In [23j, Helton 
and Vinnikov answer by the negative, showing that in the plane (n = 2) some convex 
basic semialgebraic sets cannot be LMI. An elementary example is the so-called TV screen 
set defined by the Fermat quartic 

{x G : l-xl-x^>0} (3) 

see Figure [H 
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1.1 Rigid convexity 



Assume that the set T defined in ([T]) has a non-empty interior, and choose a point Xq in 
this interior, i.e. 

Xo G int JF = {x : F{x) >- 0} 

where >- means positive definite. A segment starting from xq attains the boundary of 
when the determinant po{x) = det F{x) vanishes. The remaining polynomial inequalities 
Pi{x) > 0, i > only isolate the convex connected component containing xq. This moti- 
vated Helton and Vinnikov J23j to study semialgebraic sets defined by a single polynomial 
inequahty 

V = {xeW : p{x) > 0}. (4) 

The set {x : p{x) > 0} is called an algebraic interior with defining polynomial p{x), and 
it is equal to int V when V is convex. With these notations, the question addressed in 
[23] is as follows: what are the conditions satisfied by a polynomial p{x) so that V is an 
LMI set ? 

For notational simplicity we will assume, without loss of generality, that xq = 0, so that 
V contains the origin, and hence we can normalize p{x) so that p{0) = 1. 

If p{x) = det F{x) for some matrix mapping F{x) we say that p{x) has a determinan- 
tal representation. In particular, the polynomial po{x) in ([2]) has a symmetric linear 
determinantal representation. 

Consider an LMI set as in ([1]) and define 

p(x) = det F(x) 

as the determinant of the symmetric pencil F{x). Note that degp = m, the dimension of 
F[x). Define the algebraic variety 

C = {x G M" : p{x) = 0} (5) 

and notice that the boundary of JF is included in C. Indeed, a point x* along the boundary 
of is such that the rank of F{x*) vanishes. Since the origin belongs to JF it holds Fq >z 0. 

Now consider a line passing through the origin, parametrized as x(t, z) = tz where t G M 
is a parameter and 2 G is any vector with unit norm. For all z, the symmetric 
matrix F{x{t, z)) = Fq + t{ziFi -!-■■■ + ZnFn) has only real eigenvalues as a pencil of t, 
and its determinant t \—>- p{x{t,z)) = det F{x{t)) has only real roots. Therefore, a given 
polynomial level set P as in (j4]) is LMI only if the polynomial t i— > p{x{t,z)) has only 
real roots for all z, it must satisfy the so-called real zero condition [23]. Geometrically 
it means that a generic line passing through the origin must intersect the variety ([5]) at 
m = degp real points. The set V is then called rigidly convex, a geometric property 
implying convexity. 

A striking result of [23] is that rigid convexity is also a sufficient condition for a polynomial 
level set to be an LMI set in the plane, i.e. when n = 2. For example, it can be checked 
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easily that the TV screen set ([3]) is not rigidly convex since a generic line cuts the quartic 
curve only twice. 



In the litterature on partial differential equations, polynomials satisfying real zero condi- 
tion are also called hyperbolic polynomials, and the corresponding LMI set is called the 
hyperbolicity cone, see [36j for a survey, and ^SlJ for connections between real zero and 
hyperbolic polynomials. 

In passing, note the fundamental distinction between an LMI set (as defined above) and 
a semidefinite representable set, as defined in [331 [5j. A semidefinite representable set is 
the projection of an LMI set: 

n n„ 

J^={xEW : 3uEW- : F{x,u) = Fq + ^x^F^ + ^m^G^- h 0} 



where the variables Uj, sometimes called liftings, are instrumental to the construction of 
the set through an extended pencil F{x,u). Such a set is called a lifted LMI set. It 
is convex semialgebraic, but in general it is not basic. However, it can be expressed as 
a union of basic semialgebraic sets. In the case of the TV screen set ([3]) a lifted LMI 
representation follows from the extended pencil 



F(x, u) 



l+Ui 
U2 



U2 
1 — Ml 



i Xi 

Xi Ui 
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X2 U2 



obtained by introducing two liftings. It seems that the problem of knowing which convex 
semialgebraic sets are semidefinite representable is still mostly open, see [30i l24j for recent 
developments. 



1.2 Determinantal representation 



Once rigid convexity of a plane set, or equivalently the real zero property of its defining 
polynomial, is established, the next step is constructing an LMI representation. Alge- 
braically, given a real zero bivariate polynomial p{xi^X2) of degree m, the problem consists 
in finding symmetric matrices Fq, Fi and F2 of dimension m such that 

X2) = det(Fo + FiXi + F2X2) 

and Fq ^ 0. If the Fj are symmetric complex- valued matrices, this is a well-studied 
problem of algebraic geometry called determinantal representation, see p7] for a classical 
reference and [3], |3l] for more recent surveys and extensions to trivariate polynomials. 

If one relaxes the dimension constraint (allowing the Fj to have dimension larger than m) 
and the symmetry constraint (allowing the Fj to be non-symmetric), then results from 
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linear systems state-space realization theory (in particular linear fractional representa- 
tions, LFRs) can be invoked to design computer algorithms solving constructively the 
determinantal representation problem. For example, the LFR toolbox for Matlab [2T] is 
a user-friendly package allowing to find non-symmetric determinantal representations: 

» Ifrs xl x2 

» f=l/(l-xl"4-x2"4) 

LFR-object with 1 output(s), 1 input(s) and state(s). 

Uncertainty blocks (globally (8 x 8)): 

Name Dims Type Real/Cplx Full/Seal Bounds 

xl 4x4 LTI r s [-1,1] 

x2 4x4 LTI r s [-1,1] 

The software builds a state-space realization of order 8 of the transfer function f{x) = 
l/p{x). This indicates that a non-symmetric real pencil F{x) of dimension 8 could be 
found that satisfies detF(a;) = p{x), as evidenced by the following script using the Sym- 
bohc Math Toolbox (Matlab gateway to Maple): 

» syms xl x2 

» D=diag( [ones(l,4)*xl ones(l ,4)*x2] ) ; 
» F=eye(8)-F.a*D 
F = 



1, 


-xl. 


0, 


0, 


0, 


0, 


0, 


0] 


0, 


1, 


-xl, 


0, 


0, 


0, 


0, 


0] 


0, 


0, 


1, 


-xl, 


0, 


0, 


0, 


0] 


xl, 


0, 


0, 


1, 


-x2. 


0, 


0, 


0] 


0, 


0, 


0, 


0, 


1, 


-x2. 


0, 


0] 


0, 


0, 


0, 


0, 


0, 


1, 


-x2. 


0] 


0, 


0, 


0, 


0, 


0, 


0, 


1, 


-x2] 


xl. 


0, 


0, 


0, 


-x2, 


0, 


0, 


1] 



» det(F) 
ans = 

-x2~4+l-xl"4 

Note that LFR and state-space realization techniques are not restricted to the bivariate 
case, but they result in pencils of large dimension (typically much larger than the degree 
of the polynomial), and there is apparently no easy way to reduce the size of a pencil. 

If one insists on having the Fi symmetric, then results from non-commutative state- 
space realizations can be invoked to derive a determinantal representation, at the price of 
relaxing the sign constraint on Fq. An implementation is available in the NC Mathematica 
package Here too, these techniques may produce pencils of large dimension. 



Now if one insists on having symmetric Fi of minimal dimension m, then two essentially 
equivalent constructive procedures are known in the bivariate case to derive Hermitian 
complex- valued Fi from a defining polynomial p{xi,X2) of degree m. Real symmetric 
solutions must then be extracted from the set of complex Hermitian solutions. 

The first one is based on the construction of a basis for the Riemann-Roch space of 
complete linear systems of the algebraic plane curve C given in ([5]). The procedure is 
described in [12]: one needs to find a curve of degree m — l touching C at each intersection 
point, i.e. the gradients must match. The algorithm is illustrated in [32]. It is not clear 
however how to build a touching curve ensuring Fq y 0. 

The second determinantal representation algorithm is sketched in [23] and in much more 
detail in [15]. It is based on complex Riemann surface theory [201 IB]- Explicit expres- 
sions for the Fi matrices are given via theta functions. Numerically, the key ingredient 
is the computation of the period matrix of the algebraic curve and the corresponding 
Abel-Jacobi map. The period matrix of a curve can be computed numerically with the 
algcurves package of Maple, see [11] and the tutorial [10] for recent developments, in- 
cluding new algorithms for explicit computations of the Abel-Jacobi map. A working 
computer implementation taking p{xi,X2) as input and producing the Fi matrices as 
output is still missing however. 



1.3 Contribution 

The focus of this paper is mostly on computational methods and numerical algorithms. 
The contribution is twofold. 

First in Section [2] we describe an algorithm for detecting rigid convexity in the plane. 
Given a bivariate polynomial p{xi,X2), the algorithm uses a hybrid symbolic-numerical 
method to detect whether the connected component of the sublevel set (jlj) containing 
the origin is rigidly convex. The problem boils down to deciding whether a univariate 
polynomial matrix is positive semidefinite. This is a well-known problem in linear systems 
theory, for which numerical linear algebra algorithms are available (namely eigenvalue 
decomposition), as well as a (more expensive but more flexible) semideflnite programming 
formulation. 

Then in Section [3] we describe an algorithm for solving the determinantal representation 
problem for algebraic plane curves of genus zero. The algorithm is essentially symbolic, 
using Bezoutians, but it assumes that a rational parametrization of the curve is available. 
The idea behind the algorithm is not new, and can be traced back to [28], as surveyed 
recently in [27] . An algorithm for detecting rigid convexity of a connected component 
delimited by such curves readily follows. 

Extensions to positive genus algebraic plane curves and higher dimensional sets are men- 
tioned in Section HI In particular we survey the case of cubic plane curves and cubic 
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surfaces which are well understood. The case of quartic (and higher degree) curves seems 
to be mostly open, and computer implementations of determinantal representations are 
still missing. Similarly, checking rigid convexity in higher dimensions seems to be com- 
putational challenging since it amounts to deciding whether a multivariate polynomial 
matrix is positive semidefinite. 



2 Detecting rigid convexity in the plane 

In this section we design an algorithm to assess whether the connected component delim- 
ited by a bivariate polynomial around the origin is rigidly convex. The idea is elementary 
and consists in formulating algebraically the geometric condition of rigid convexity of the 
set V defined in a line passing through the origin cuts the algebraic curve C defined 
in ([S]) a number of times which is equal to the total degree m of the defining bivariate 
polynomial 

X e p{x) = ^ PaX" = Poo + PioXi + P01X2 + P20XI + PnXlX2 H 

aeN"^ ,\a\<m 

A line passing through the origin can be parametrized as: 

Xi = r cos 9 = {z~^ + z) . . 

X2 = r sin 9 = it~^{z~^ — z) 

where z = e*^ and t = 2r~^. Along this line, we define 

m 

t G M ^ q{t) = rp{x) = J2 (lk{z)t'' (7) 

k=0 

as a univariate polynomial of degree m which vanishes on C. Moreover q{t) is monic since 
<im{z) = p{0) = 1. The remaining coefficients are Laurent polynomials 

m 
fc=0 

with real coefficients, also called trigonometric cosine polynomials. Set V is rigidly convex 
if and only if this polynomial has only real roots, i.e. if the number of intersections of the 
line with the curve C is maximal. 



2.1 Counting the real roots of a polynomial 

A well-known result of real algebraic geometry [21 Theorem 4.57] states that a univariate 
polynomial q{t) of degree m has only real roots if and only if its Hermite matrix is positive 
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semidefinite. The Hermite matrix is the m-hy-m moment matrix of a discrete measure 
supported with unit weights on the roots xi, . . . ,Xm of polynomial q{t) (note that these 
roots are not necessarily distinct). It is a symmetric Hankel matrix whose entries are 
Newton sums A^i+j = ^^Li ^I^"'- The Newton sums are elementary symmetric functions of 
the roots that can be expressed explicitly as polynomial functions of the coefficients of q{t). 
Recursive expressions are available to compute the Nk, or equivalently, Nk = trace C'^ 
where C is a companion matrix of polynomial q{t), i.e. a matrix with eigenvalues Xi, 
see e.g. [21 Proposition 4.54]. Recall that coefficients of the polynomial q{t) given in 
(I7j) are Laurent polynomials. It follows that the Hermite matrix of q{t) is a symmetric 
trigonometric polynomial matrix of dimension m, that we denote by H{z). We have 
proved the following result. 

Lemma 1 The bivariate polynomial p{x) is rigidly convex if and only if its Hermite 
matrix H{z) is positive semidefinite along the unit circle. Coefficients of H{z) are explicit 
polynomial expressions of the coefficients ofp{x). 



2.2 Positive semidefiniteness of polynomial matrices 

The problem of checking positive semidefiniteness of a polynomial matrix on the unit 
circle is generally referred to as (discrete-time) spectral factorization. It is a well-known 
problem of systems and circuit theory [ITJ HHl [19]. The positivity condition can also 
be defined on the imaginary axis (continuous-time spectral factorization) or the real axis. 
Various numerical methods are available to solve this problem [2H]. Several algorithms are 
implemented in the Polynomial Toolbox for Matlab [35]. In increasing order of complexity, 
we can distinguish between 

• Newton-Raphson algorithms: the spectral factorization problem is formulated as a 
quadratic polynomial matrix equation which is then solved iteratively [26] . At each 
step, a linear polynomial matrix equation must be solved ^5j. Quadratic (resp. 
linear) convergence is ensured locally if the polynomial matrix is positive definite 
(resp. semidefinite); 

• polynomial operations: a sequence of elementary operations is carried out in the 
ring of polynomials to reduce the polynomial matrix to some canonical form, see [8] 
and [48j for a recent survey. These algorithms are cheap computationally but their 
numerical behavior (stability) is unclear; 

• algebraic Riccati equation: using state-space realization, the problem is formulated 
as a quadratic matrix equation, which in turn can be solved via a matrix eigenvalue 
decomposition with a particular structure [Ml Uni 1^ : 

• semidefinite programming: polynomial matrix positivity is formulated as a con- 
vex semidefinite program, see [H] and the recent surveys [171 H] • The particular 
structure of this semidefinite program can be exploited in interior-point schemes. 
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in particular when forming the gradient and Hessian. General purpose semidefinite 
solvers can be used as well. 



The semidefinite programming formulation of discrete-time polynomial matrix factoriza- 
tion, a straightfoward transposition of the continuous-time case studied in [H], is as 



follows. The symmetric trigonometric polynomial matrix H(z) = Hq + Hi{z + z ) + 
■ ■ ■ + Hd{z'^ + z~'^) of size m is positive semidefinite along the unit circle if and only if 
there is a symmetric matrix P of size dm such that 





' Ho 


Hi ■ 


■ Hd' 






L{P) = 


Hi 








+ 




/ 




Hd 


■ 


■■ 




■ 


■ 



p 



' I 


" 




I 



" Ho 


Hoi 








+ 







P[B\A] 



(8) 



P[D\C ]hO. 





I 


T 


I 


B^{z'^)L{P)B{z) = 


z-^ 


L{P) 


z 











Notice that the columns and rows of the above matrix are indexed w.r.t. increasing powers 
of z in such a way that 



Hiz) 



Positive semidefiniteness of L{P) then amounts to the existence of a polynomial sum- 
of-squares decomposition of H{z). From the Schur decomposition L{P) = U'^U with 
U = [Uo Ui ■■■ Ud] it follows that 

H{z) = U{z-YU{z). (9) 

Polynomial matrix U{z) = Uo + Uiz + ■ ■ ■ + Udz'^ is called a spectral factor. 

If the LMI problem ([H]) is feasible, then it admits a whole family of solutions. Assuming 
that Ho >~ ^1 maximizing the trace of P subject to the LMI constraint ([8]) yields a 
particular solution P* such that rank L{P*) = m. It follows that the Schur complement 
of 

L{P) 



Ho 


+ B^PB - 


D^PD 


-k 


. Hqi 


+ A^PB - 


C^PD 


A^PA - C^PC 



w.r.t. Hq vanishes, where symmetric entries are denoted by -k. This means that P* satisfies 
the quadratic matrix equation 

A^PA - C^PC - {H^^ + A^PB - C^PD)Hq\Hoi + B'^PA - D'^PC) = 

called the (discrete-time) algebraic Riccati equation. In this case, the spectral factor U (z) 
in ([9]) is square non-singular. 
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2.3 Example: cubic curve 



Consider the component of set (j4]) around the origin dehmited by the cubic polynomial 
p{x) = 1 — Xi — Ax\ — X2 + 4xf, see Figure [2l 




Figure 2: Cubic curve and its component around the origin (shaded). 



Using the substitution ([6]) we obtain 

q{t) = 12{z + z-^) + A{z^ + z-^) - (10 + 3{z^ + z-^))t - {z + z-^)t^ + . 
From the companion matrix 



C 



z + z-'^ 10 + 3(^2^^-2) _i2(z + 2-1) -4(^3^ 
1 

1 



we build (symbolically) the Hermite matrix 
3 



Hiz) 



z + z 



2 I ^-2-] 



22 + 7{z^ + z 



•k 

-k 



22 + 7(^2 + ^-2) 6(z + 2^1) -2(^3 + ^-3) 250 + 124(^2 ^ ^-2) _^ ^5(^4 ^_ ^-4^ 



in 



Solving (numerically) the LMI ([8]) with SeDuMi interfaced with YALMIP yields the spec- 
tral factorization ([9]) with factor (in Matlab notation) 



Uiz) 



-0.9021 -0.7094^2 -0.5284^ + 0.20272^ -11.7639 - 9.6359^^ - 1.5201^^ 

0.1925^ 4.3449 + 1.6218^2 0.7771z - 0.5411;z3 

1.1578 -0.5527^2 0.3819;z + 0.1579^^ 2.4331 - 2.8689^^ - 1.1844^^ 



which certifies numerically that p{x) is rigidly convex. 



2.4 Example: quartic curve 

Let us apply the algorithm to test rigid convexity of the TV quartic level set ([3]) with 
p{x) = 1 — xf — X2, see Figure [1] 

We obtain q{t) = —12 — 2{z^ + z~^) + 1'^ and the Hermite matrix 

4 * -k -k 

★ * 

^^^ 48 + 8(^^ + ^-4) k ■ 

_ 48 + 8{z^ + z-^) _ 

From the zero diagonal entries and the non-zero entries in the corresponding rows and 
columns we conclude that H{z) cannot be positive semidefinite and hence that the TV 
quartic level set is not rigidly convex. 



2.5 Numerical considerations 



Since the Hermite matrix H{z) has a Hankel structure, and positive definite symmetric 
Hankel matrices have a conditioning (ratio of extreme eigenvalues) which can be bounded 
below by an exponential function of the matrix size [H US] , it may be appropriate to apply 
a congruence transformation on matrix H{z), also called scaling. 

For example, if H{e^^°) is positive definite for some ^0 (say Oq = 0, but other choices are 
also possible), it admits a Schur factorization H{e^^°) = V^DV with V orthogonal and D 
diagonal non-singular. If D is reasonably well-conditioned, we can test positive semidefi- 
niteness of the modified trigonometric polynomial matrix Hq{z) = V D~^/'^H{z)D~^^'^V^ 
along the unit circle, which is such that iJo(e*^°) is the identity matrix. If D is not well- 
conditioned, we can still use Ho{z) = VH{z)V^ which is such that ifo(e*^°) is a diagonal 
matrix. 

The impact of this data scaling on the numerical behavior of the semidefinite programming 
or algebraic Riccati equation solvers is however out of the scope of this paper. 
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3 LMI sets and rational algebraic plane curves 



In the case that the algebraic curve C in ([5]) has genus zero, i.e. the curve is rationally 
parametrizable, an alternative algorithm can be devised to test rigid convexity of a con- 
nected component delimited by C. The algorithm is based on elimination theory. It 
uses a particular symmetric form of a resultant called the Bezoutian. As a by-product, 
the algorithm also solves the determinantal representation problem in this case. As sur- 
veyed recently in [27], the key idea of using Bezoutians in the context of determinantal 
representations can be traced back to ^28j . 

Starting from the implicit representation 

C = {x eR^ : pix) = 0} (10) 

of curve C, with p{x) a bivariate polynomial of degree m, we apply a parametrization 
algorithm to obtain an explicit representation 

C = {x e : xi = qi{u)/qo{u), X2 = q2{u)/qo{u), m G R} (11) 

with qi{u) univariate polynomials of degree m. Algorithms for parametrizing an implicit 
algebraic curve are described in [U [3H1 US]. An implementation by Mark van Hoeij is 
available in the algcurves package of Maple. The coefficients of qi{u) are generally found 
in an algebraic extension of small degree over the field of coefficients of p{x). 

With the help of resultants, we can eliminate the variable u in parametrization ( fTTj) and 
recover an implicit equation ( fTOj) . see P Section 3.3]. To address this implicitization 
problem, we make use of a particular resultant, the Bezoutian, see [15], Section 5.1.2]. 
Given two univariate polynomials g,h of the same degree m (if the degree is not the 
same, the smallest degree polynomial is considered as a degree m polynomial with zero 
leading coefficients) build the following bivariate polynomial 

u — V ^ ^ 

fc=0 z=o 

called the Bezoutian of g and /i, and the corresponding symmetric matrix B{g,h) of size 
m X m with entries bki bilinear in coefficients of g and h. As shown e.g. in [ini Section 
5.1.2], the determinant of the Bezoutian matrix is the resultant, so we can use it to derive 
the implicit equation (fTOj) of a curve from the explicit equations ffTTj) . 



Lemma 2 Given polynomials qo, qi, q2 in 177]) . a polynomial p in / fiOj) is given by p{x) = 
detF(x) where 

F{x) = B{qi,q2) + XiB{q2,qo) + X2B{qi,qo) . . 

= F0 + X1F1 + X2F2 ^ ' 

is a symmetric pencil of size m. 
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Proof: Rewrite the system of equations ( fTTl) as 



92{.u) 



q2{u) -X2qo{u) 








and use the Bezoutian resultant to ehminate indeterminate u and obtain conditions for a 
point (xi, X2) to belong to the curve. The Bezoutian matrix is B{gi, §2) = B{qi —Xiqo, q2 — 
X2qo) = B{qi, ^2) + XiB{q2, go) + 2^2-8(^1, qo)- Linearity in x follows from bilinearity of the 
Bezoutian and the common factor go 

Lemma [2] provides an implicit equation of curve (fTOl) in symmetric linear determinantal 
form. 



3.1 Detecting rigid convexity 

Once polynomial p{x) is in symmetric linear determinantal form as in Lemma [2], checking 
rigid convexity of the connected component containing the origin x = amounts to testing 
positive definiteness of F{0) = Fq = B{qi, q2). 

Lemma 3 The Bezoutian matrix B{qi,q2) is positive semidefinite if and only if polyno- 
mial qi{u) and q2{u) have only real roots that interlace. 

Proof: The signature (number of positive eigenvalues minus number of negative eigen- 
values) of the Bezoutian of qiiu) and q2{u) is the Cauchy index of the rational function 
qi{u)/q2{u), the number of jumps of the function from —00 to +00 minus the number of 
jumps from +00 to —00, see |2l Definition 2.53] or [21 Theorem 9.4]. It is maximum when 
B{qi, q2) is positive definite. This occurs if and only if the roots of qiiu) and q2{u) are all 
real and interlace. □ 

Lemma 4 The connected component around the origin delimited by curve l[W\) is rigidly 
convex if and only if B{qi, q2) >z 0. 

Proof: Since Fq >z 0, the set admits the LMI representation {x G : F{x) y 0}, which 
is equivalent to being rigidly convex. □ 



3.2 Finding a rigidly convex component 

If the connected component around the origin is not ridigly convex, it may happen that 
there is another rigidly convex connected component elsewhere. To find it, it suffices to 
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determine a point x* such that F{x*) >z 0. This is equivalent to solving a bivariate LMI 
problem. 

We can apply primal-dual interior-point methods [33] to solve this semidefinite program- 
ming problem, since the function f{x) = — \ogp{x) = log det F(a;)~^ is a strictly convex 
self-concordant barrier for the interior of the LMI set. If the LMI set is bounded, min- 
imizing f{x) yields the analytic center of the set. If the LMI set is empty, the dual 
semidefinite problem yields a Farkas certificate of infeasibility. However, in the bivariate 
case a point x* satisfying F[x*) >z can be found more easily with real algebraic geometry 
and univariate polynomial root extraction. 

A first approach consists in identifying the local minimizers of function f{x). They are 
such that the gradient g{x) of p{x) vanishes, i.e. they are such that 

g^(^:,) = ^ = tmce{p{x)F-\x)F,) = 0. (13) 

OXi 

We can characterize these minimizers by eliminating one variable, say xi, from the system 
gi{x) = g2{x) = 0, and solving for the other variable X2 via polynomial root extraction. 
Resultants can be used for that purpose. 

A second approach consists in finding points on the boundary of the LMI set, which are 
such that p{x) = and either gi{x) = or g2{x) = 0. Here too, resultants can be applied 
to end up with a polynomial root extraction problem. 

From the points generated by these two procedures, we keep only those satisfying F{x) >z 
0, an inequality that can be certified by testing the signs of the coefficients of the charac- 
teristic polynomial of F{ explained in the introduction. 



3.3 Example: Capricorn curve 



Let p{x) = xl{x1 + X2) — 2(xf + X2 — X2y. With the parametrization function of the 
algcurves package of Maple, we obtain a rational parametrization 

= 45 _ 8t + lOt^ + 

qi{t) = -7 + Ut - 18^2 - + t4 
g2(t) = 49 _ 28t - lOt^ + 4^3 + 

With the BezoutMatrix function of the Linear Algebra package, we build the correspond- 
ing symmetric pencil 



F{x) 



8 - 4x1 - 4x2 
8 + 20x1 - 28x2 
-72 + 20x1 + 52x2 
56 - 4x1 - 52x2 



40 + 60x1 + 92x2 
-8 - 36x1 - 84x2 776 4 
-168 + 180x1 + 180x2 -952 



540x1 + 476x2 
- 940x1 + 740x2 



1960 - 868x1 - 1924x2 
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whose determinant (up to a constant factor) is equal to p{x). The eigenvalues of F{0) are 
equal to (double) and 1392 ± 48v533. They are all non-negative which indicates that 
the origin lies on the boundary of an LMI region defined by F{x) >z 0. 



Values of X2 at local optima satisfying the system of cubic equations (fT3|) can be found 
with Maple as follows: 

> p:=xl~2*(xl"2+x2"2)-2*(xl"2+x2"2-x2)"2: 

> solve (resultant (dif f (p,xl) ,dif f (p,x2) ,xl) ) ; 

0, 0, 0, 1, 1/2, 3+sqrt(5), 3-sqrt(5), 3+sqrt(5), 3-sqrt(5) 

from which it follows that, say, the point xi = 0, X2 = 1/2 is such that F{x) >- 0. 




Figure 3: Capricorn curve defining an LMI region (shaded). 



The corresponding LMI region together with the quartic Capricorn curve are represented 
on Figure [31 



3.4 Example: bean curve 



0.8 



0.6 - 



0.4 - 



0.2 - 



-0.2 - 



-0.4 - 



-0.6 - 




-0.8' ' ' ' ' ' ' ' 

0.2 0.4 0.6 0.8 1 

Figure 4: Bean curve defining a region which is not LMI (shaded). 
Let p{x) = xf + x\x2 + X2 — x\+xix1. With the Bezoutians we obtain the following pencil 



F{x) = 



Xi 




-k 


-k 


X2 


1 


-k 


-k 


Xi 







-k 


X2 


1 — Xi 





1 - 



We can check that -F(O) has eigenvalues 2 and (triple) and this is the only point for 
which F{x) ^0. It follows that the convex set delimited by the curve p{x) = is not 
LMI, see Figure H 



4 Extensions 

In this paragraph we outline some potential extensions of the results to algebraic plane 
curves of positive genus and varieties of higher dimensions. 



4.1 Cubic plane curves 



The case of cubic plane algebraic curves is well understood, see e.g. [H] or [30]. Singular 
cubics (genus zero) can be handled via Bezoutians as in Section [31 Smooth cubics (genus 
one), also called elliptic curves, can be handled via their Hessians. 

Let p{xi, X2) be a cubic polynomial that we homogeneize top{xo, xi, X2) = x'qp{xi/xo, X2/X0) 
Define its 3-by-3 symmetric Hessian matrix H{p{x)) with entries 

^ ^ d'^pjx) 
dxidxj 

and the corresponding Hessian h{x) = det H{p{x)). The elliptic curve p{x) = has 9 in- 
flection points, or flexes, satisfying p(a;) = h{x) = 0, and 3 of them are real. Since p(x) and 
h{x) share the same flexes and the Hessian matrix yields a symmetric linear determinantal 
representation for h{x), we can use homotopy to flnd a determinantal representation for 
p{x). 

For real t deflne the parametrized Hessian g{x,t) = det H{h{x) + tp{x)) and flnd t* 
satisfying g{x*,t*) = p{x*) at a real flex x* by solving a cubic equation. As a result, 
we obtain three distinct symmetric pencils not equivalent by congruence transformation. 
One of the may be deflnite hence LMI. 



Build the Hessian h{x) = det H{p{x)) = 8(xq + 

2 

XqX ^ 



For example, let p{x) = x^ — x^ — Xi. 

SxqxI — 3x1X2) and the parametrized Hessian g{x,t) = det H{h{x) +tp{x)) 
57Qt'^xlxi + ■■■ + 110592a;f. Polynomial g{x,t) matches g{x) at flex Xq 
{0, 24, —24} yielding the following three representations 



24^3^ 
for t* G 



[X 



1 

Xi 



-Xi -k 

1 



F'^ix) 



F^{x) 



4-3 



4-^ 



1 + 3x1 






-X2 


—1 — Xi 


-k 


-1 + Xi 


-X2 


1 — Xi 


1 — 3xi 


-k 


-k 


-X2 


1 — Xi 


-k 


1 + Xi 


X2 1 


+ Xi 



such that detF*(x) = p{x) for all i = 1,2,3. Only the flrst one generates an LMI set 
FHx) y 0. 
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4.2 Positive genus plane curves 



The case of algebraic plane curves of positive genus and degree equal to four (quartic) 
or higher is mostly open. Whereas rigid convexity of higher degree polynomials can be 
checked with the proposed approach, there is no known implementation of an algorithm 
that produces symmetric linear determinantal (and hence LMI) representations in this 
case. For quartics, contact curves can be recovered from bitangents. In [13] complex 
symmetric linear determinantal representations of the quartic 1 + + x| could be derived 
from the equations of the bitangents found previously by Cayley for this particular curve. 



Bezoutians can be generalized to the multivariate surveyed in [15]. In Lemma [2] 

we derived a symmetric linear determinantal representation by eliminating the variable u 
in the system of equations 

gi{u) = qi{u) - Xiqo{u) = 
g2{u) = q2{u) - X2qo{u) = 

corresponding to a rational parametrization Xi{u) = qi{u)/qo{u), X2{u) = q2{u)/qo{u) of 
the curve p{xi, X2) = 0. In the positive genus case, such a rational parametrization is not 
available, but we can still define a system of equations 

gi{ui,U2) = xi-ui =0 

92{Ul,U2) = X2-U2 =0 
93iUi,U2) = p{Ui,U2) = 

after eliminating variables ui and U2- Define the discrete 



describing the curve X2) 
differentials 

g{ui,U2) - g{vi,V2) 



dig{u,v) 

and the quadratic form 
gi digi d2gi 



Ui - Vi 



d2g{u,v) 



g{vi,U2) - g{vi,V2) 



U2 - V2 



det 



92 
93 



di92 d2g2 
digs d2gz 



det 



Xi — Ui 
X2 - U2 
P{Ul,U2) 



-1 



dip{u, v) 





-1 
d2p{u, v) 



a,/3 



using bi-indices a and (3. Then the matrix F{x) of the quadratic form is a symmetric pencil 
satisfying det F{x) = p{x)q{x) where q{x) is an extraneous factor. We hope that q{x) does 
not depend on x, even though this cannot be guaranteed in general. For example, in the 



case of the Fermat curve p{x) = 1 — x^ 
package for Maple [6], we could obtain 



whose genus is three, using the mult ires 



F{x) 



which is such that det F fx) 



-1 














Xi 


X2 











Xi 





-1 

















X2 





-1 





Xi 





-1 

















X2 





-1 








Xi 


-1 

















X2 





-1 














p{x) 


, i.e. 


q{x) 




1. 
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4.3 Surfaces and hypersurfaces 



The case n = m = 3, i.e. cubic surfaces, is well understood, see [7] for a full constructive 
development. All the self-adjoint linear determinantal representations can be obtained 
from the tritangent planes. The number of non-equivalent representations depends on 
the number and class of real lines among the 27 complex lines of the surface. See [39] for 
a nice survey on cubic surfaces. 




Figure 5: The Cayley cubic surface with its convex connected component. 



A well-known example is the Cayley cubic 

1111 
— + — + — + — = 

Uo Ui U2 U3 



1Q 



whose algebraic equation is 

UQU1U2 + U0U1U3 + U0U2U3 + U1U2U3 = 0. 
Under involutary linear mapping 







■ 1 


1 


1 


1 




Uo 


Xi 


1 


1 


1 


-1 


-1 




Ui 


X2 


~ 2 


1 


-1 


1 


-1 




U2 


X3 




1 


-1 


-1 


1 




U3 



the dehomogenized {xq = 1) algebraic equation becomes 





1 


Xi 


X2 




p{x) = 1 — xl — xl — xl — 2x1X2X3 = det 


Xi 


1 




= det F{x 




X2 


X3 


1 





which is the determinant of the 3x3 moment matrix of the MAXCUT LMI relaxation. 
The surface p{x) = is represented on Figure [5l using the surf visualization package. In 
particular, we can easily identify the convex connected component containing the origin, 
described by the LMI F{x) >z 0. The component has four vertices, or singularities, for 
which the rank of F{x) drops down to one. 

In general, only curves and cubic surfaces admit generically a determinantal representa- 
tion. When n > 3 or m > 3 and no lifting is allowed, the hypersurface p{x) = must be 
highly singular to have a determinantal representation [3], and hence, a fortiori, an LMI 
representation. This leaves however open the existence of alternative algorithms consisting 
in constructing symmetric linear determinantal representations of modified polynomials 
p{x)q{x), with q{x) globally nonnegative, say q{x) = for > 1 large 

enough. 

Finally, let us conclude by remarking that, as a by-product of the proof leading to Lemma 
[H checking numerically rigid convexity of a scalar polynomial when n > 2 amounts to 
checking positivity of a multivariate Hermite matrix. See e.g. [13] for recent developments 
on the use of semidefinite programming for multivariate trigonometric polynomial matrix 
positivity. 
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